A Semigroup Theory of Rate Independent Hysteresis 
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We explore a macroscopic, algebraic approach to rate independent hysteresis using semigroup 
theory. A macroscopic description of metastable states relevant to rate independent hysteresis is 
introduced using field history. The semigroup structure of the history space is identified. Using 
semigroup theory and related mathematical techniques, the general relation between return point 
memory (RPM) and partial order is discovered. For hysteresis system with RPM, a variational 
principle is identified. The erasing properties of field histories are also characterized. The connection 
between this semigroup approach and other models are also discussed. 

PACS numbers: 75.60.Ej, 03.65.Fd, 02.90.+p 
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Rate independent hysteresis refers to irreversible phe- 
nomena in the adiabatic limit. It is universally exhib- 
ited by disordered, frustrated systems at low tempera- 
ture, such as random magnets, glasses, plastic materials, 
etc. These systems typically have exponentially large 
number of metastable states. When thermal fluctuations 
are negligible, such a system stay indefinitely in whatever 
metastable state it happens to find. Consequently, the 
system properties depend on the whole history of exter- 
nal field back to infinite past. 

Theories of rate independent hysteresis [l|, [3, H, |3| 
has been studied traditionally by the engineer commu- 
nity. Engineers have long understood that bi-stability (or 
rather multi-stability) is the essential ingredient of hys- 
teresis, and have proposed a wide range of models [22] to 
account for these phenomena, such as the play model, the 
stop model, the Preisach model the Stoner-Wohlfarth 
model 6]. These models are usually based on vari- 
ous simple microscopic mechanisms and provide definite 
schemes to calculate hysteresis loops. Mathematicians 
have also came tip with models involving partial differ- 
ential equations [3| . 

The physics of zero temperature disordered system has 
also attracted considerable interests in physics commu- 
nity in recent years. One primary interests seems to 
be the zero temperature critical phenomena exhibited 
by hysteresis loops of random field Ising model (RFIM) 
0, Q and other related models , as one tunes dis- 

order variance continuously. Another class of system is 
depinning of elastic manifolds inside a random potential 
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12l |. The physics of the pinned phase, however, has 
not been systematically explored, largely due to the lack 
of appropriate formalism and language. 

Nevertheless, the general macroscopic theory of rate 
independent hysteresis is largely undeveloped. In the au- 
thor's opinion, such a theory should care less about mi- 
croscopic mechanisms responsible for hysteresis and focus 
on the macroscopic, material independent aspects: How 
does a hysteretic system remember its history? What is 
the macroscopic meaning of "memory" ? How can mem- 
ory be erased? Exploration of these issues not only is im- 



portant to advance our understanding of nonequilibrium 
physics, it may also lead to new designs and applications 
in information science, material science and engineering. 
In this letter, we introduce such a macroscopic formalism 
involving semigroup theory, as well as a few exact results 
about rate independent hysteresis and return point mem- 
ory. Detailed discussions and proofs will be presented in 
a separated publication [l§| • 

We shall only consider hysteresis systems driven by 
a scalar external field, for the sake of simplicity. To 
develop a macroscopic approach to hysteresis, we first 
need a macroscopic description, or a "name" , for each 
metastable state relevant to rate independent hysteresis. 
For a system at thermal equilibrium, a state is fully char- 
acterized by a few macroscopic parameter, such as vol- 
ume, entropy, or their conjugate external field, such as 
pressure, temperature, etc. For a disordered system at 
zero temperature, these parameters do not uniquely de- 
termines the current state. This is of course the essence 
of hysteresis: the system remembers its past. The first 
important observation is that for a given hysteresis dy- 
namics, and a given initial state, a history of external field 
corresponds to a unique final state. Therefore we can la- 
bel a metastable state by the field history that we use to 
prepare this particular state. That is, the field history 
is the "order parameter" of hysteresis systems. Two im- 
portant comments: 1) There are metastable states that 
are not reachable by any field history. However, these 
states are not relevant to the study of rate independent 
hysteresis. 2) More importantly, there may be multiple 
field histories that correspond to the same state. This is 
the case of hysteresis system with RPM. We will show 
that the redundancy in description can be "gauged out" 
inside the framework of semigroup theory. Study of this 
"gauge symmetry" help us to understand the basic al- 
gebraic structure of hysteresis phenomena with return 
point memory. 

Let us start from some initial metastable state at a 
field value ft-o- Very often we choose the initial state to 
be one of two saturated states |-f) or |— ), at the positive 
or negative saturation fields. Km or —Hm respectively. 
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FIG. 1: A: The external field (solid line) as a function of time 
determines a metastable state, for a given initial state. Ar- 
bitrary time reparameterizations of this function represents 
the same field history. The dashed line is another field his- 
tory which is "smaller" than the first one, according to the 
relation > defined in the history space. B: New mechanism 
for information storage and read-out using rate independent 
hysteresis. 

[2^ We may introduce two type of operators, r]{Ah), 
which increase the field by Ah, and ^(A/i), which de- 
crease the field by Ah. A metastable states reached by 
a particular hysteresis process then can be represented 
as a string of operators 77 and ^ acting on the initial 
state. For example, for the history shown in Fig. [TJ^, 
we start from the negatively saturated state |— ) at field 
—hM and raise, lower, and raise the field, end up with a 
state with final field value /13. This state can be repre- 
sented as \4>) =Tl{h-i + h2)i{h2+hi)r){hi+hM)\-). [2I 
The string of operators ri{h^ + /i2)»7(/i2 + hi)ri(hi + hu) 
then represents the particular field history we used to 
prepare the metastable state 1-0). Note that 1) raising or 
lowering the field by 5h always changes the metastable 
states, no matter how small 5h is; 2) The field history and 
the initial state determines at which field the final state 
is metastable. For example, the state shown in left panel 
of Fig. [l]is metastable at field ~hM + (^1 + h^i) ~ + 
hi) + (/13 -1-/12) = We can naturally put two histo- 
ries together and form a product. The product of ry(A/ii) 
and ^(A/12) is just ^(A/i2)?7(Aft,i), which amounts to first 
increase field by A/ii, then decrease it by A/i2- This is 
not the same as 7]{Ahi — A/12), as the dynamics is irre- 
versible. It is also different from r/(A/ii)^(A/i2), exactly 
for the same reason. Therefore the history product is not 
commutative. On the other hand, the product of i]{hi) 
and rj{h2) is just rj{hi + ^12), which follows from adia- 
baticity. Product of more complicated field histories can 
be recursively defined. This product is clearly associa- 
tive, i.e. A{BC) = {AB)C. Therefore the set of all field 
histories forms a semigroup [l5| . which we shall call 
the general history semigroup. Now take a state Al-^o} 
obtained by applying a field history A onto the initial 
state IV'o)- Suppose we want to do hysteresis experiment 
by applying operation B onto the system, the new state 
of the system after this experiment is simply given by 
BA\iIiq). That is, the rate independent hysteresis dynam- 
ics is simply represented as the product (concatenation) of 
histories. This strongly suggests that semigroup theory 
is the natural language for rate independent hysteresis. 

Almost all of the hysteresis models studied up to now 



possess a very important property, that is, return point 
memory (RPM, or wipe out property): If one raises the 
external field by h, then lowers it by h, and finally raises 
it again by h, then end result on the system is identi- 
cal to that if one raises the field by h and does nothing 
else. RPM implies the closeness of hysteresis subloops 
regardless of the initial state. Sethna and collaborators 
|16| discovered a deep result about RPM: If there is a 
partial order in the configuration space and the hystere- 
sis dynamics respects this order, i.e. no passing, then the 
system must exhibit return point memory. Qualitatively 
speaking, a partial order is a binary relation > in the 
configuration space that allows one to compare certain 
pairs of state without inconsistency. It has to satisfy: 1) 
reflexivity: a > a; 2) antisymmetry: if a > b and b > a, 
then a = b; 3) Transtivity: if a > 6 and b > c then a > c. 
A dynamics in the configuration space satisfies the no 
passing rule, if it respects the partial order. That is, if 
we have two initial states \ipi{Q)) > |'02(O)), which are 
driven by two field histories hi{t) > h2{t) respectively. 
Then for arbitrary time t > 0, we have |'0i(i)) > |V'2(0)- 
Qualitatively, RPM means that a system starting behind 
and driven by a lower field can not pass another system 
starting in front and driven by a higher field. This no- 
passing property was first explored by Middleton 17[ , in 
the setting of driven charge density waves. It was sub- 
sequently realized that it applies to a wide class of zero 
temperature driven systems. 

The relation between partial order and return point 
memory discovered by Sethna et al is elegant, deep, and 
general. It illustrates how general properties of hystere- 
sis follow from a few essential features of the dynamics. 
One naturally wonders whether the converse is true: If 
RPM holds, does it imply a partial order in the config- 
urational space that is respected by the dynamics? As 
a special case, recently Narayan et al found that ID 
antifcrromagnetic Ising chain shows RPM if starts from 
saturation state, but not if starts from a generic state. 
Even though the dynamics does not respect the trivial 
partial order defined by spin orientation, Narayan found 
a hidden partial order that satisfies no passing rule. A 
general answer to this question can be obtained from the 
semigroup approach introduced in this work. 

Return point memory is represented by the following 
identities in the general history semigroup: 

m^mih) = vmwvih) = ^lih). (i) 

It defines an equivalence relation in the general history 
semigroup. All histories equivalent to each other under 
Eq. ID) correspond to the same metastable state. The 
history semigroup for hysteresis system with RPM is then 
the quotient semigroup over the equivalence Eq. ([1]). The 
description of metastable state in terms of field history 
therefore contains redundancy: return point memory is 
a gauge symmetry in the history space. 
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We can define a binary relation > in the history space: 
let A and B be two field histories, which are strings of 
operators 77's and ^'s. We define A > B (" A is larger 
than B") if, when starting from the same initial field, 
there exists a time parameterization of these field his- 
tories such that A{t) > B{t) for all t. An example is 
shown in Fig. [T]A_. The relation > is clearly reflective and 
transitive. More importantly, one can prove that if the 
return point memory Eg. (fT|) holds, then the relation > 
is also antisymmetric [l3|, i.e. it is a partial order. Fur- 
thermore, one trivially see this partial order is respected 
by the hysteresis dynamics, i.e. no passing rule holds: If 
A > B, and C > D, then AC > BD by definition of 
the relation > . Combining this result with that of refer- 
ence [3l , we arrive at the following basic theorem of rate 
independent hysteresis: 

Theorem (RPM and partial order): a rate inde- 
pendent hysteresis system exhibits return point memory 
if and only if there exists a partial order in the space of 
all reachable metastable states that is respected by the 
dynamics (i.e. no-passing). 

We emphasize that the partial order > is defined only 
in the history space, which is in one-to-one correspon- 
dence to the set of metastable states that are reachable 
by rate independent hysteresis processes. This is a sub- 
space of the whole configurational space. [2^ It may not 
be related to the naive partial order (defined by spin ori- 
entations in spin models for example) in any apparent 
fashion. Neither we should not expect that this partial 
order can be extended to the whole configurational space. 
One recent example is given by Narayan et al [l8| . 

The one-to-one correspondence between field histories 
and metastable states reachable by a hysteresis process 
may provide a novel mechanism of information storage. 
We can imagine to encode some data into the extremal 
values of a field history in a rate independent fashion, 
which is in turn applied to a zero temperature hysteretic 
system, such as a disordered magnet. This process im- 
prints the data sequence onto the hysteresis system. One 
important question is then how to extract informations 
stored in the system. For a system with return point 
memory, there is a simple mechanism achieving this goal: 
One applies a field which is monotonically increasing over 
time. Whenever the field passes an extreme value of 
the original field history, one should observe a slope dis- 
continuity in the input-output curve. This mechanism, 
schematically illustrated in Fig. [T}3, allows one to re- 
cover all the local maximums of the original field history. 
Further research along this direction may lead to feasi- 
ble applications of hysteresis properties in information 
science and technology. 

There are other problems that we need to worry about. 
The information stored in a hysteresis system may get 
erased accidently by noise. On the other hand, if we 
want to imprint new information onto a system, we want 
to make sure that the old information can be completely 



(-hM, hM) 



P hm 



(-hw , h| 



-hi, hi) 



1 


m 















a 


P ' 


-1 




A 





(h,,h,) a 



FIG. 2: A: Input-output relation for a hysteron. h is the exter- 
nal field, m is the "magnetization". B: Staircase representa- 
tion of Preisach model and any hysteresis system with RPM. 
The partial order of field histories is also represented graphi- 
cally as the relation "above" and "below" between staircases. 

erased. Again, these questions can be properly addressed 
using semigroup theory. We define another binary rela- 
tion l> [erasing) in the history space: A history A erases 
another history B, or At> B,ii and only if AB = A. We 
see that the field history B has no effect if it is followed 
by the history A, that is, B is erased by A. The following 
two results follow from basic properties of semigroup: 

• All histories that erase the history A form a sub- 
semigroup, which we call the eraser subsemigroup 
of A. 

• All histories that are erased by the history A form 
a subsemigroup, which we shall call the erased sub- 
semigroup of A. 

The structures of eraser semigroup and erased semigroup 
of a given history should be problems of fundamental 
importance in the general theory of hysteresis, since they 
charcterize the macroscopic meaning of memory. For a 
general hysteresis system without RPM, one can show 
that the only way to erase a history is to raise or lower 
the field to saturation. On the other hand, sweeping the 
field to saturation clearly erase all field histories. For 
system with RPM, the author managed to completely 
characterize the structures of the eraser semigroup and 
the erased semigroup of a special subclass of field history. 
Due to limit of space, results are not presented here but 
will be published elsewhere [isj . 

The Preisach model ^ is one of the most successful 
phenomenological models for rate independent hysteresis. 
It models hysteresis system as a superposition of many 
independent hysterons, i.e. basic units of hysteresis with 
two switch fields a < /3, as illustrated in Fig. [2K. A hys- 
teron exhibits bi-stability when a < h < (3. It flips up 
(its value changes from -1 to 1) when the external field 
passes f3 from below, and fiips down (its value changes 
from 1 to -1) when the field passes a from above. The 
Preisach model admits a partial order (defined by the ori- 
entation of all hysterons) and no-passing rule, from which 
the return point memory follows. Mayergoyz [20| proved 
that if a hysteresis system satisfies return point memory 
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and congruence property, then it can be represented by 
the Preisach model. 

The Preisach model admits an extremely useful 
graphic representation as shown in Fig. [1)3 [2l|: In the 
plane of two switch fields a and every point inside 
the triangle represents a hysteron. A metastable state 
reachable by hysteresis processes starting from satura- 
tion states is then represented by a staircase running 
from one of the catheti to the hypotenuse. All hysteron 
above/below the staircase are turned down/up. For a 
given distribution of hysterons, this graphic representa- 
tion constitutes a convenient starting point for the cal- 
culation of hysteresis loops [H in the input-output plane. 
This staircase representation also gives an intuitive mean- 
ing to the partial order: a metastable state A is larger 
than another one B if and only if the staircase of A is 
one the top of B everywhere. Finally the hysteresis dy- 
namics of raising/lowering the field is represented by the 
deformation of the staircase to the above/left [2]] |. 

We make another observation that the staircase is a 
faithful representation of field history for any hysteresis 
system with RPM. Using this graphic representation, one 
can show the following variational principle of rate 
independent hysteresis with RPM: Let \4>) be an 
arbitrary initial state that is metastable at field value 
ho. Then the state 7]{h)\ip) is the lower bound of all 
states that are 1) metastable at field value ho + h, and 2) 
larger than {ip). Similarly the state £,{h)\ip) is the upper 
bound of all states that are 1) metastable at field value 
ho — h, and 2) smaller than lip). This variational princi- 
ple completely defines the hysteresis dynamics, just like 
the Hamiltonian variational principle completely defines 
the classical Lagrange dynamics. It is amazing that an 
irreversible process such as rate independent hysteresis 
can be formulated from a variational principle. 

What is the difference between the Preisach model 
and the semigroup theory of hysteresis with RPM? The 
semigroup theory addresses the most general properties 
of rate independent hysteresis that are independent of 
microscopic details. The semigroup theory itself does 
not give prediction of the shapes of hysteresis loops. 
For the latter purpose we need some constitutive rela- 
tions. In the Preisach model, these are given by the 
independent-hysteron assumption and the hysteron dis- 
tribution. Naturally different materials are characterized 
by different constitutive relations. In a more general ver- 
sion of Preisach model, it is perceivable to give up the 
independent hysteron hypothesis and consider (weak) in- 
teractions. This has indeed been studied by Mayergoyz 
and others [H, 0]. It should be emphasized that disen- 
tanglement of the material-independent hysteresis prop- 
erties from the materials-dependent ones is an important 
step towards a better understanding of hysteresis theory, 
and this can be achieved by the semigroup approach. 
Also further study of the variational principle may lead 
to discovery of more systematic approximation methods 



for the calculation of material properties. It may also 
help us identify the connection between lattice models 
such as RFIM and more phenomenological models such 
as Preisach models. These will be explored in the future. 
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